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基于 因 于 分 析 的 简 文 染 桥 损 伤 识 别 


字 , 张 素 侠 , 刘 习 军 
(天 津 大 学 机 械 工程 学 院 ,300072 天 津 ) 


要 :为 了 降低 测量 噪声 对 桥梁 损伤 检测 结果 的 影响 ,提高 检测 方法 在 实际 工程 中 的 实用 性 , 提 


出 了 一 种 基于 因子 分 析 的 简 支 梁 桥 损伤 识别 方法 。 首 先 , 对 采集 的 位 移 信号 进行 因子 分 析 , 得 到 信 
号 的 公共 因子 和 特殊 因子 ,通过 去 除 公共 因子 达 目的 。 然 后 ,采用 统计 方法 和 D-S 证 据 理 论 
对 特殊 因子 进行 处 理 并 构造 损伤 指标 。 通 过 一 个 梁 桥 有 限 元 模型 的 算 例 进行 数值 分 析 , 讨 论 


了 嗓 声 、 车 辆 移动 速度 及 质量 对 损 pr o 0 车 速 在 30 ~80 km/h 范围 内 时 ， 


“所 提 方 法 可 以 准确 定位 较 小 程度 的 结构 损伤 ;此 外 , 骂 声 和 载 质 量 对 识别 结果 影响 较 小 ,满足 实际 


工 况 需 求 ,说 明 所 提 方 法 具有 很 强 的 抗 品 能 力 和 实用 价值 。 
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Damage identification of simply supported beam 
bridge based on factor analysis 


WANG Zhiyu,ZHANG Suxia, LIU Xijun 
(School of Mechanical Engineering, Tianjin University,300072 Tianjin, China) 


3Abstract: In order to reduce the influence of measurement noise on the detection results of bridge damage 


and improve the practicability of the detection method in practical engineering,a damage identification 


method of simply supported beam bridge based on factor analysis was proposed. Firstly ,factor analyses of 


the collected displacement signals were made to obtain the common factors and special factors of the sig- 


nals ,and achieved the purpose of noise reduction by removing the common factors. Then ,statistical meth- 


ods and D-S evidence theory were used to process the special factors and construct damage indicators. The 


influences of noise ,vehicle moving speed and vehicle mass on damage identification were explored through 


numerical analysis of a finite element model of a simply-supported beam bridge. The results showed that 


when the vehicle speed was in the range of 30-80 km/h,the method could accurately locate the structural 


damage with a small degree. In addition ,the effect of noise and vehicle mass on the identification results 


was small , which met the requirements of the actual working conditions ,indicating that the method had a 


strong anti-noise ability and practical value. 
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近年 发 展 起 来 的 基于 结构 振动 的 动力 学 损伤 检 
测 方法 , 因 其 具有 操作 简便 、 经 济 实用 等 优点 ,在 结 
构 损伤 检测 领域 得 到 了 广泛 的 应 用 。 马 爱 敏 等 建 
立 了 裂纹 简 支 梁 连 续 抗 弯 模 型 ,并 通过 求解 裂纹 梁 
的 动力 指纹 ,包括 固有 频率 、 振 型 以 及 振 型 曲率 ,说 
明了 基于 动力 指纹 检测 裂纹 损伤 的 可 行 性 。 唐 盛 华 
等 5 利用 结构 损伤 前 后 的 柔 度 矩阵 曲率 范 数 差 成 功 
识别 了 连续 梁 桥 的 损伤 。 李 雪 艳 等 ” 将 模式 匹配 法 
与 结构 的 位 移 响应 协 方差 参数 变化 比值 相 结合 ,对 
结构 的 损伤 位 置 进行 了 识别 。 该 方法 无 需 结构 模 
型 ,计算 比较 简便 。 

然而 在 实际 工程 中 ,结构 的 动态 响应 必然 会 受 
到 噪声 ,温度 等 外 界 干扰 信息 的 影响 ,甚至 掩盖 原本 
的 损伤 信息 。 外 界 环境 因素 的 影响 已 经 限制 了 
很 夫 一 部 分 损伤 识别 方法 在 具体 工程 中 的 实现 。 
本 为 了 减少 噪声 的 影响 ,一 些 学 者 采用 扩展 卡尔 
曼 肖 波 算法 对 信号 进行 处 理 。Sen 等 中 针对 卡尔 晤 
滤 纺 算法 在 特定 情况 下 受到 噪声 影响 较 大 的 问题 ， 
提名 采用 双 扩 展 卡 尔 曼 滤波 技术 对 结构 参数 和 状态 
进 每 同步 估计 ,并 通过 数值 模型 证 明了 该 方法 的 适 
甩 俐 。 何 浩 翔 等 中 将 扩展 卡尔 曼 滤波 与 静 力 凝 至 法 
结合 起 来 ,利用 静 力 凝聚 法 等 效 简化 复杂 梁 桥 的 模 
型 $ 然 后 采用 扩展 卡尔 曼 滤 波 成 功 识 别 了 三 跨 连 续 
梁 焰 的 损伤 。Lai 等 中 提出 移动 窗口 的 扩展 卡尔 曼 
滤波 算法 。 相 对 传统 算法 而 言 ,该 方法 对 初始 值 不 
壬 感 ,在 损伤 检测 方面 具有 更 强 的 鲁 棒 性 , 对 系统 参 
数 鲍 识别 也 更 准确 。 虽 然 目前 一 些 关 于 扩展 卡尔 曼 
滤波 的 研究 可 以 在 降低 噪声 干扰 的 影响 下 成 功 识别 
结 情 的 损伤 ,但 是 由 于 方程 中 状态 向 量 维 数 较 大 以 
及 结构 状态 和 参数 向 量 之 间 的 非 线性 看 全 ,计算 比 
较 繁琐 , 且 容 易 出 现 不 适 定 问题 ;所 以 要 运用 在 实际 
工程 中 还 需要 进一步 研究 。 

为 了 能 够 在 未 知 噪 声 及 其 他 环境 因素 的 情况 下 


分 析 等 非 参 数 化 方法 对 数据 进行 处 理 。 常 鹏 等 ” 利 
用 主 成 分 分 析 法 减少 了 温度 对 损伤 特征 参数 的 影 
响 ,并 结合 小 波 包 变 换 构 造 损 伤 指标 对 结构 进行 损 
伤 识别 。Wah 等 "认为 结构 振动 特征 参数 的 改变 
是 结构 损伤 和 环境 因素 的 共同 作用 ,并 成 功 采用 主 
成 分 分 析 法 减 小 了 环境 因素 的 影响 。 刁 延 松 等 
提出 了 基于 AR 模型 和 因子 分 析 的 结构 损伤 识别 方 
法 。 首 先 利用 因子 分 析 法 降低 温度 对 AR 模型 系数 
的 影响 ,然后 结合 标准 差 控 制图 对 结构 进行 损伤 识 
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别 。 然 而 ,上 述 考 虑 环境 因素 的 研究 工作 仅仅 利用 
一 些 方法 降低 了 环境 因素 对 结构 动力 响应 特征 参数 
的 影响 ,对 后 续 的 损伤 识别 方法 并 未 做 过 多 的 研究 ， 
比较 依赖 于 现 有 的 方法 ,导致 整个 识别 过 程 比 较 
繁琐 。 

为 了 在 降低 噪声 影响 的 同时 提高 检测 方法 的 实 
用 性 ,本 研究 提出 了 一 种 基于 因子 分 析 的 简 支 梁 桥 
损伤 识别 方法 。 首 先 ,采用 因子 分 析 法 将 结构 的 位 
移 响 应 分 为 公共 因子 和 特殊 因子 两 部 分 ,通过 去 除 
公共 因子 来 降低 测量 噪声 的 影响 ;然后 ,结合 D-S 证 
据 理论 ,通过 融合 计算 特殊 因子 的 统计 特征 参数 来 
构造 损伤 指标 。 数 值 算 例 中 考虑 了 实际 高 速 路 况 下 
的 车 速 范 围 以 及 工程 测量 时 的 噪声 强度 大 小 。 算 例 
分 析 和 模型 实验 的 结果 说 明 本 研究 所 提 方 法 不 仅 可 
以 准确 定位 桥梁 结构 的 损伤 , 而 且 在 抵抗 测量 噪声 
方面 具有 很 大 的 优势 ,有 很 强 的 实用 价值 。 


1 因子 分 析 


因子 分 析 首 先 需 要 计算 出 各 变量 的 均值 和 标准 
差 来 对 变量 进行 标准 化 计算 ,然后 通过 计算 标准 化 
后 变量 的 相关 系数 矩阵 ,将 每 一 个 样本 的 多 个 复杂 
变量 表示 成 少数 几 个 公共 因子 的 线性 组 合 与 特殊 因 
子 之 和 。 其 中 ,公共 因子 代表 同一 样本 所 有 变量 
含 的 共同 信息 ,而 特殊 因子 则 表示 每 一 个 变量 所 独 
有 的 信息 。 

对 于 一 个 B 维 数据 处 ,假设 其 具有 7/ 个 公共 因 
子 , 则 通过 因子 分 析 计 算 , 可 以 将 表示 为 


Xl CI Qa Qn A 21 
2 Ci2 022 az | 万 22 

有著 = = . -+| -|(1) 
Xp Qlp U2B Wp f EB 


式 中 :六 ,六 为 数据 于 的 /个 公共 因子 ;ax 为 第 
了 个 公共 因子 在 站 的 第 j 个 分 量 % 上 的 载荷 ;s; 
为 % 的 特殊 因子 。 

式 (1) 也 称 为 因子 模型 ,可 以 将 其 写 为 


X=Af+e (2) 
其 中 
Ql 0Q21 Qn 
a a a 
12 22 2 
A= 
Up QBp “** Qip 
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= [六 fi] ,se = [si ,ea 26] 
式 中 ,4 为 因子 载荷 矩阵 。 将 式 (2) 写 为 特殊 因子 的 
表达 式 , 即 

=Af-X (3) 

从 式 (3) 可 以 看 出 ,只 要 计算 出 原始 数据 中 每 

个 分 量 的 特殊 因子 , 即 可 减少 原始 数据 中 公共 因子 

的 成 分 。 而 损伤 检测 过 程 中 的 噪声 温度 等 环境 因 

素 会 对 结构 的 整体 产生 影响 ,存在 于 数据 的 每 个 变 

量 中 ,属于 公共 因子 部 分 ,所 以 通过 提取 特殊 因子 可 
以 达到 减少 环境 因素 影响 的 目的 下。 


2 损伤 指标 的 构建 


一 采集 移动 荷载 作用 下 桥梁 结构 个 测 点 的 位 移 
响应 如 下 。 


和 yi yy2 全 y? 
© Ea (4) 
yy ye 


式 噶 :n, 为 单个 测 点 的 采样 点 数 ;y 为 结构 在 第 i 个 
测 列 采集 的 第 人 个 时 刻 的 位 移 响应 。 
中 本 研究 分 两 部 分 进行 指标 的 计算 。 
= 1) 以 测 点 为 样本 ,时 间 为 变量 的 因子 分 析 计算 。 
-过 针对 式 (4) ,以 测 点 作为 样本 ,时 间作 为 变量 ， 
可 入 得 到 ”个 样本 及 刀 个 变量 ,对 这 种 情况 下 的 数 
据 和 进行 因子 分 析 计算 。 为 了 满足 因子 分 析 的 前 提 检 
验 么 件 (KMO 和 Barttett 球 度 检验 ) , 可 以 采取 删除 
某 此 相关 性 较 强 的 变量 的 方法 进行 计算 ,最 后 提取 
出 特殊 因子 e, 。 该 部 分 的 特殊 因子 表示 同一 测 点 在 
不 同时 刻下 响应 的 差异 ,其 统计 特征 参数 具有 凸显 
损伤 信息 的 能 力 。 

第 i 个 测 点 的 特殊 因子 的 方差 Vs; ,和 峰 度 系数 
Kss,; 分 别 为 

Vsp ; = 二 [a = Dy 2 (5) 


ti k=1 


1 
| 
ee Elik— ~ 
7 £21 ' 


Ky; = A 

gp,i a, (6) 
式 中 ,e114 表示 结构 第 i 个 测 点 第 个 时 刻 位 移 响 应 
的 特殊 因子 。 


当 梁 上 各 个 相 邻 传感器 之 间 的 距离 相等 时 ,分 
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的 方差 变化 率 Ssrv 和 峰 度 系数 变化 率 Ssrx 的 绝对 值 。 


Vs i — Vs i 

[Sr | = ee (7) 
Le i — 人; 4 

[Sani ee (8 ) 


式 中 ,h 为 两 相 邻 测 点 的 间距 。 为 了 综合 Swv 和 Sspx 
的 绝对 值 进行 损伤 识别 ,需要 用 D-S 证 据 理论 对 其 
进行 融合 。 

D-S 证 据 理 论 通 过 计算 每 个 证 据 的 基本 概率 分 
配 函 数 ( BPA) ,并 结合 Dempster 合成 规则 来 达到 毕 
合 考 虑 多 个 证 据 并 给 出 最 终结 果 的 目的 。 其 中 ， 
BPA 也 称 为 mass 函数 ,用 m 表示 。 在 损伤 识别 方 
面 ,可 设 所 有 识别 结果 组 成 的 集合 为 9, 即 辨识 框 
架 。 若 4 为 @ 的 子 集 且 相应 的 基本 概率 分 配 函 数 
m(4) 关 0 时 , 则 称 4 为 m 的 焦 元 。 现 有 辨识 框架 下 
的 2 个 证 据 分 别 为 和 ,对 应 的 mass 函数 分 别 
为 mm 和 7m ,对 应 的 焦 元 分 别 为 4 和 4;, 根 据 合 成 规 


则 ,可 得 融合 后 的 证 据 集合 C 为 
0 C= 
m(C) = De (As) ma A 人 
一 CO,C 
外 = oy mi(A;)m,(4,) 


(9) 
要 融合 Sssy 和 Ssrx ,需要 计算 二 者 的 mass 函数 
mi 和 m, ,把 结构 的 每 个 测 点 c; 作为 D-S 证 据 理 论 
的 子 集 , 则 计算 式 为 
| | 


> | SS 
i=1 


oe 
m,(c,) = | Sree 


2 Ssec.s| 
(10) 
结合 式 (9) 和 式 (10) 对 Ssv 和 Ssx 的 绝对 值 进 
行 融 合计 算 ,将 计算 结果 记 为 D,。 
2 ) 以 时 间 为 样本 , 测 点 为 变量 的 因子 分 析 计 算 。 
对 于 式 (4) , 以 时 间作 为 样本 , 测 点 作为 变量 ， 
可 得 到 ,个 样本 及 nn 个 变量 。 对 数据 进行 因子 分 
析 计 算 并 提取 相应 的 特殊 因子 e,。 这 部 分 的 特殊 因 
子 的 含义 为 同一 时 刻 结构 在 不 同 测 点 下 响应 的 差 
异 ,其 本 身 就 合 有 结构 的 损伤 信息 。 
第 i 个 测 点 的 特殊 因子 的 偏 度 系 数 Ss, 可 以 表 


mi(c:) = 


别 按照 式 (7) 和 式 (8) 计 算 相 邻 2 个 测 点 间 特 殊 因子 


示 为 


2 ( 至 > Bh 

对 > Cw 2 ex | 

各 个 测 点 的 特殊 因子 的 中 位 数 用 SFM 表示 。 
根据 式 (7) 和 式 (8) 的 方法 ,可 以 计算 出 相 邻 2 个 测 
点 间 特 殊 因子 的 偏 度 系数 变化 率 Ssps 以 及 中 位 数 变 
化 率 Ssmv 的 绝对 值 。 

同样 地 ,结合 式 (9) 对 Ssss 和 5sm 的 绝对 值 进行 
D-S 融合 计算 ,二 者 mass 函数 的 算法 与 式 (10) 相 
同 ,将 计算 结果 记 为 D,。 

定义 最 终 用 于 损伤 检测 的 损伤 指标 D; 为 数值 
D, 和 D, 进 行 D-S 证 据 融 合 后 的 数值 的 绝对 值 。 


Ssp,i 二 (11) 


3 > 数值 算 例 分 析 


3 了 算 例 模型 
家 采用 如 图 1 所 示 的 简 支 梁 桥 有 限 元 模型 进行 数 


广 上 ,~V, 人 人 入 甘 日 工 
CNINAAIVD 1 FA 
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值 动力 分 析 计 算 。 桥 全 长 24 m, 全 桥 划 分 为 30 个 单 
元 ,并 在 桥 的 下 方 均匀 布置 9 个 位 移 传感器 。 定 义 1 
号 梁 段 为 图 1 中 从 左 向 右 第 1 个 传感器 和 第 2 个 传 
感 器 之 间 的 部 分 ( 即 单元 4 ~ 单元 6) ,其 他 梁 段 号 以 
此 类 推 ,已 在 图 1 中 标 出 。 相 关 参 数 设置 如 下 :弹性 
模 量 =3.45 x10” Pa, 桥 梁 横 截面 积 为 1.5 m , 截 
面 惯性 矩 1=0. 281 25 m ,密度 p =2500 kg/m , 泊 松 
比 A =0.2。 以 弹簧 质量 块 模 型 来 模拟 车 辆 荷载 , 质 
量 块 重 10 kN, 弹 得 的 刚度 系数 为 1 x10” N/m, 和 车辆 
的 行驶 速度 v=10 m/s 。 假 定 结 构 的 损伤 只 与 单 
元 刚度 变化 有 关 , 并 通过 单元 刚度 折 减 来 模拟 损伤 。 
共 设 置 3 种 损伤 工 况 ,具体 如 表 1 所 示 。 
表 1 数值 模拟 的 损伤 工 况 


Tab.1 Damage cases of numerical simulation 


损伤 工 况 损伤 单元 号 损伤 程度 /% 
1 16 5 
2 24 5 
3 14,24 53 


0.8mX30=24m 


人 位 移 传感器 损伤 单元 


简 支 梁 模 型 及 加 载 方式 示意 图 


Fig.1 Schematic of the simply-supported bridge and loading mode 


3725 损伤 识别 结果 


对 采集 到 的 位 移 信 号 分 别 以 时 间 和 测 点 作为 样 
本 进行 因子 分 析 计 算 ,得 到 2 种 情况 下 的 特殊 因子 ， 
并 采用 本 研究 第 2 章 的 指标 构造 方法 进行 计算 ,得 
到 损伤 识别 指标 D;。 本 研究 因子 分 析 计 算 均 通过 
SPSS 软件 进行 。 经 过 计算 ,本 研究 数值 模型 输出 数 
据 的 KMO 检验 值 均 在 0.8 以 上 ,适合 做 因子 分 析 。 
此 外 , 公 因 子 数目 的 确定 可 根据 特征 值 大 于 1 以 及 
方差 累积 解释 大 于 85% 的 原则 ,最 终 确定 的 公 因 子 
个 数 为 2。 

为 了 更 好 地 表示 出 损伤 位 置 ,采用 设置 预警 线 
的 方式 ,计算 方法 如 下 


(12) 


= 

V 
三 | 一。 
i 
局 


Sax 二 人 ? Sin 去 1 ( 13 ) 


式 中 :M; 为 第 i 段 梁 的 损伤 指标 值 ;MM 为 各 梁 段 损 
伤 指标 的 最 大 值 ;N 为 梁 段 总 数 ;$ 为 安全 系数 。 当 
式 (12) 成 立时 ,认为 第 工段 梁 发 生 损伤 。 安 全 系数 
5 越 小 判断 条 件 越 保 守 , 其 取 值 应 该 介 于 最 大 安全 
系数 Su 和 最 小 安全 系数 Su 之 间 。 二 者 可 以 按照 
式 (13 ) 进行 计算 ,本 研究 中 的 $ 取 二 者 的 中 间 值 。 

按 上 述 方法 计算 出 工 况 1 和 工 况 2 下 损伤 指标 
的 预警 线 , 并 绘 出 相应 的 损伤 识别 结果 图 ,如 图 2 所 
示 。 这 2 个 工 况 均 为 单 损伤 工 况 , 图 2(a) 中 5 号 梁 
的 指标 值 较 大 ,对 应 于 16 单元 损伤 ,图 2(b) 中 7 号 
梁 的 指标 值 较 大 ,对 应 于 24 单元 的 损伤 。 识 别 结 晤 
与 预 设 工 况 相符 。 
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开 (b) 工 况 2 
这 。 图 2 工 况 1 和 工 况 2 的 损伤 识别 结果 


CmFig. 2 Damage identification results of cases 1 and 2 
中 为 了 更 清晰 地 描述 本 研究 方法 的 识别 过 程 ,下 
工 况 3 为 例 进行 说 明 。 图 3 为 跨 中 5 号 测 点 的 
便 咎 响应 ,分 别 以 时 间 和 测 点 为 样本 ,对 所 有 测 点 的 
啊 霹 数据 进行 因子 分 析 计算 ,可 得 各 测 点 标准 化 后 
的 全 因子 重 构 位 移 响应 。 用 标准 化 后 的 位 移 响应 减 
共 兴 因子 重 构 响应 , 即 可 得 各 测 点 的 特殊 因子 。 然 
后 禄 据 第 2 章 的 方法 将 不 同情 况 下 特殊 因子 的 统计 
循 铂 参数 进行 融合 计算 ,计算 结果 如 图 4 所 示 。 其 
中 四 ,为 以 测 点 为 样本 计算 得 到 的 值 。 可 以 看 出 ， 
根 纪 D, 只 能 识别 7 号 梁 的 损伤 ，D, 为 以 时 间 为 样 
本 和 鸯 计算 结果 ， 虽 然 图 4 (b) 中 4 号 梁 和 7 号 梁 
处 的 D, 值 都 比较 大 ,但 是 考虑 到 预警 线 的 值 ， 也 
只 能 判定 4 号 梁 一 处 损伤 ; 图 4 (c) 为 D, 和 D, 融 
合 后 的 识别 结果 ， 图 中 损伤 梁 段 处 的 指标 值 D; 均 
明显 超过 了 预警 线 ， 所 以 能 够 根据 记 准确 判 断 出 
结构 的 损伤 位 置 。 
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友 3 测 点 5 的 位 移 响 应 ( 工 况 3) 


Fig.3 Displacement response of measuring point 9 (case 3) 
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人 六 


二 二 损伤 指标 


梁 段 号 
(0) D, 
图 4 Di、D,、D; 的 损伤 识别 结果 ( 工 况 3) 


Fig.4 Damage identification result of Di ,D, and Di (case 3) 


3.3 ”噪声 对 损伤 识别 的 影响 

为 了 验证 本 研究 所 提 方 法 的 降 品 效果 ,在 结构 
的 位 移 响应 中 添加 高 斯 白 噪 声 ,其 计算 公式 可 表 
未 为 


yi,noise = Yi + E,Nioisestd( y;) (14) 
式 中 :yi ,wise 为 添加 噪声 后 的 动态 位 移 响应 ;Ek, 为 品 
声 水 平 ;NN,, 为 均值 为 0, 标 准 差 为 1 的 随机 数 序列 ; 


std(y;) 为 位 移 响应 的 标准 差 。 此 外 ,这 里 添加 的 品 
声 是 用 来 模拟 测量 噪声 的 , 而 不 只 是 环境 噪声 ,还 包 
括 路 面 粗糙 度 .不 平顺 等 其 他 环境 因素 对 响应 造成 
的 干扰 。 

以 工 况 2 为 例 ,分 别 在 响应 中 添加 噪声 水 平 为 
10% .20% .30% 的 高 斯 白 噪声 ,损伤 识别 结果 如 图 5 
所 示 。 可 以 看 出 , 随 着 噪声 强度 的 逐渐 增 大 ,未 损伤 
梁 段 的 损伤 指标 值 也 有 所 增加 ,但 变化 幅度 较 小 且 
始终 没有 超过 所 设 定 的 预警 线 , 而 损伤 梁 段 (7 号 
梁 ) 的 指标 值 一 直 远 超 预 警 线 ,说 明 本 研究 所 提出 的 
损伤 识别 方法 具有 很 好 的 噪声 鲁 棒 性 。 
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CN 图 5 噪声 影响 下 的 损伤 识别 结 
CN Fig.5 Damage identification results 
~ under the influence of noise 


3 小 ”荷载 移动 速度 对 损伤 识别 的 影响 


开 趟 同 荷 裁 速度 的 作用 下 , 桥 面 所 受 的 激励 不 同 ， 
所 采集 到 的 动力 响应 信号 也 不 同 。 为 了 研究 荷载 移 
动 媳 度 对 损伤 识别 的 影响 ,以 工 况 2 为 例 , 将 荷载 的 
移动 速度 分 别 设 为 8 m/s 和 12 m/s, 并 计算 出 相应 
的 损伤 指标 D;。 计 算 结 果 如 图 6 所 示 。 

对 比 图 2(b) 与 图 6 可 以 看 出 ,虽然 所 提出 的 损 
伤 识别 方法 在 不 同 的 荷载 移动 速度 下 均 能 准确 识别 
损伤 位 置 ;但 与 图 2(b) 和 图 6(b) 相 比 ,图 6(a) 中 未 
损伤 梁 段 的 损伤 指标 更 加 接近 预警 线 ,对 损伤 识别 
的 干扰 程度 较 大 。 

为 了 进一步 研究 荷载 速度 对 损伤 识别 的 影响 ， 
定义 除 损伤 位 置 外 的 其 他 节点 处 损伤 指标 最 大 的 点 
为 干扰 点 ,并 定义 干扰 点 和 损伤 处 的 损伤 指标 的 比 
值 为 干扰 率 ,干扰 率 越 大 ,说 明 识别 效果 越 不 理想 。 
按照 这 种 方法 ,计算 出 上 述 计算 结果 的 干扰 率 。 此 
外 ,为 了 探究 干扰 率 与 荷载 移动 速度 之 间 的 关系 ,并 
且 考 虑 实际 高 速 路 桥 的 车 速 范围 ,增加 计算 荷载 速 
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度 为 5 m/s9 ms 11 ms 17 m/s、22 m/s、27 m/s 和 
33 m/s 这 几 种 情况 下 的 干扰 率 。 最 后 ,为 了 研究 噪 
声 在 这 种 情况 下 的 影响 ,对 每 一 个 速度 下 的 位 移 响 
应 添加 噪声 水 平 为 30% 的 高 斯 白 噪 声 ,并 计算 出 相 
应 的 干扰 率 。 为 了 便于 观察 规律 ,将 上 述 计 算 结果 
绘 成 曲线 图 ,如 图 7 所 示 。 


25| ” 叶 e 损 伤 指标 D， 
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1.0 
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(b) w=12 m/s 
图 6 不 同 荷载 移动 速度 下 的 损伤 识别 结果 


Fig.6 Damage identification results 


with different moving load speeds 
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一 一 无 噪声 


0.1 最 佳 点 (10, 0.15) 
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荷载 移动 速度 /ms ) 


图 7 ”和 荷载 移动 速度 干扰 率 曲线 图 
Fig.7 Interference rate curves of moving load speed 
从 图 7 可 以 看 出 ,无 噪声 的 原始 曲线 中 ,干扰 率 

最 小 的 点 (最 佳 点 ) 对 应 的 荷载 速度 值 为 10 m/s ,此 

时 干扰 率 为 0.15。 相 对 于 该 速度 ,速度 过 大 或 过 小 均 

会 导致 干扰 率 的 增加 ,影响 损伤 识别 效果 甚至 导致 损 

伤 误 判 。 添 加 噪声 后 的 干扰 率 曲 线 具 有 相同 的 规 

律 ,这 也 验证 了 本 研究 方法 具有 良好 的 噪声 鲁 棒 性 。 
为 了 说 明 规律 的 普 适 性 ,对 原 桥 梁 有 限 元 模型 

参数 ( 即 桥梁 横 截 面积 .材料 密度 、 弹 性 模 量 等 ) 进 

行 修改 ,修改 后 的 参数 为 :弹性 模 量 3.55 x 10”Pa， 

密度 2 600 kg/m , 泊 松 比 0.2 ,截面 面积 为 1.2 m”， 

截面 惯性 矩 为 0. 144 mm 。 在 图 7 结果 的 基础 上 , 增 


738 
加 计算 荷载 速度 为 Sm/s .6 m/s .7 ms .9 ms 17 m/s、 
22 m/s ,27 m/s 和 33 m/s 情况 下 的 干扰 率 。 计 算 结 
果 如 图 8 所 示 。 


最 佳 点 (8, 0.13) 


0 10 20 30 40 
荷载 移动 速度 /(m *，s ) 


到 8 改变 模型 参数 后 的 速度 干扰 率 曲线 图 


Fig.8 Interference rate curve of load speed 


after changing model parameters 

对 比 图 7 和 图 8 可 以 看 出 , 蓓 载 移动 速度 的 干 
扰 率 确实 存在 上 述 总 结 的 规律 ,只 是 最 佳 点 的 位 置 
3 发 生变 化 。 为 了 进一步 总 结 采 用 本 人 研究 方法 准确 
识 测 结构 损伤 的 速度 范围 ,再 次 改变 桥梁 有 限 元 模 
型 殉 结 构 参 数 ,采用 相同 方法 进行 干扰 率 的 计算 。 
后 ,综合 这 3 种 情况 ,给 出 进行 损伤 识别 实验 的 建 

读 鹃 度 范围 为 in 


3 所 荷载 质量 对 损伤 识别 的 影响 


忆 褒 载 质量 不 同 会 引起 桥 染 所 受 激励 的 变化 。 为 
了 银 究 移动 荷载 质量 对 损伤 识别 的 影响 ,以 工 况 2 
为 例 , 在 3.2 节 研 究 的 基础 上 ,增加 计算 车 辆 荷载 为 
镶 t.7t 和 9+ 情 况 下 的 损伤 指标 ,并 绘 出 相应 的 识 
别 严 果 图 。 计 算 中 车 辆 移动 速度 均 设 为 10 m/s。 最 
终 的 计算 结果 如 图 9 所 示 。 

〇 从 图 9 可 以 看 出 ,虽然 随 着 荷载 质量 的 增加 ,一 
些 未 损伤 位 置 的 指标 值 有 所 增 大 ,但 并 没有 引起 误 
判 , 说 明 荷 载 质量 对 本 研究 提出 的 损伤 识别 方法 影 
响 较 小 。 为 了 进一步 研究 荷载 质量 的 影响 规律 , 同 
样 采用 计算 干扰 率 的 方法 绘 出 荷载 影响 下 的 干扰 率 
曲线 图 ,并 在 干扰 率 较 小 处 适当 加 密 计算 ,最 终结 
如 图 10 所 示 。 

从 图 10 可 以 看 出 ,荷载 质量 对 损伤 识别 的 影响 
规律 同 荷载 速度 类 似 , 即 存在 一 个 最 佳 载重 值 使 得 
干扰 率 最 小 ,数值 过 大 或 过 小 均 会 导致 识别 效果 变 
差 。 不 同 之 处 在 于 ,即使 在 荷载 质量 达到 较 大 载重 
(9 0 的 情况 下 ,干扰 率 的 值 也 小 于 0.4 ,并 不 会 导致 
损伤 误 判 。 为 了 验证 这 一 规律 的 普 适 性 ,同样 采用 
改变 桥梁 模型 参数 的 方法 ,得 到 改变 参数 后 荷载 质 
量 的 干扰 率 曲线 如 图 11 所 示 。 从 前 文 识别 结果 可 
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以 看 出 ,干扰 率 在 0.5 及 其 以 下 均 不 会 导致 损伤 误 
判 , 所 以 综合 图 10 和 图 11 可 以 说 明 , 对 于 不 同 的 桥 
梁 ,荷载 质量 对 损伤 识别 的 影响 较 小 , 且 存 在 上 述 总 
结 的 规律 。 
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Fig.9 Damage identification results with different load mass 
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Fig. 10 Interference rate curve of load mass 
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Fig. 11 


Interference rate curve of load mass after 


changing model parameters 
SP D0 人 0 人 2 
4 实验 验证 


一 为 验证 本 研究 方法 的 有 效 性 ,利用 简 支 梁 桥 模 
型 渤 行 实验 。 模 型 材料 为 有 机 玻璃 (PMMA ) , 密度 
六 5201 kg/m’ ,弹性 模 量 为 3. 75 x 10" Pal7] 。 实 验 
模型 尺寸 为 1. 40 m x0.20 mx0.01m (长 x 宽 x 
让 5 全 桥 共 划 分 为 14 个 单元 ,单元 长 度 为 0. 1 m， 
桥梁 下 方 均匀 布置 6 个 测 点 。 桥 梁 模 型 及 测 点 布置 
简 败 如 图 12 所 示 。 深 段 号 的 定义 与 第 3 章 相同 ,已 


MA12 113 114 
回 “ | 


© 
0.1 mxX14=1.4 m 
~> 激光 位 移 传感器 ”四 损伤 单元 
S< 图 12 实验 桥梁 模型 与 测 点 布置 示意 图 
& Fig.12 Schematic of experimental bridge 
SS 


model and measuring points 

斧 实 验 桥梁 模型 为 环形 桥 的 一 部 分 ,是 一 端 固定 

一 代 简 支 的 简 支 梁 桥 模型 。 沿 模型 长 度 方向 均匀 布 

置 6 个 测 点 ,每 个 测 点 下 方 放置 一 个 激光 位 移 传 感 

器 ,用 于 结构 位 移 啊 应 的 采集 ,所 用 的 信号 采集 设备 

是 DHDAS 动态 信号 采集 分 析 系 统 。 具 体 的 实验 布 
置 图 如 图 13 所 示 。 


罗 13 实验 布置 示意 图 


Fig. 13 Schematic of experimental set-up 
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对 于 实验 采用 的 矩形 截面 梁 , 刚度 的 折 减 可 以 
等 效 为 桥 面 宽度 的 减 小 ,所 以 实验 中 采用 治 桥 面 宽 
度 方向 作对 称 切口 来 模拟 损伤 ,本 次 实验 采用 宽度 
=100 mm 的 切口 ,如 图 14 所 示 。 损 伤 程度 根据 截 
面 惯性 矩 的 变化 来 计算 ,而 截面 惯性 矩 的 变化 可 以 
由 这 2 个 对 称 切口 的 深度 之 和 24d 与 桥 面 宽度 的 比 
值 进行 计算 。 共 设置 2 种 损伤 工 况 ,如 表 2 所 示 。 


EE 


14 ”损伤 单元 示意 


Fig.14 Schematic of damage unit 
表 2 实验 损伤 工 况 


Tab.2 Damage cases of experiments 


损伤 工 况 损伤 单元 号 切口 深度 wmm ”损伤 程度 /% 


1 5 10 10 


2 5,11 10,10 10,10 

实验 采用 质量 为 0.5 kg 的 小 车 模型 对 桥梁 进行 
激励 ,小 车 行驶 速度 为 0.4 m/s ,采样 频率 为 200 Hz， 
采样 时 间 为 4s。 对 2 种 工 况 下 采集 到 的 位 移 响应 
采用 本 研究 所 提 方 法 进行 处 理 。 经 过 计算 发 现 , 实 
验 数 据 的 KMO 值 也 达到 0.8 以 上 ,适合 做 因子 分 
析 。 同 样 根据 特征 值 大 于 1 以 及 方差 累积 解释 大 于 
85% 的 原则 ,确定 出 公 因子 的 数目 为 2。 最 终 的 损 
伤 识别 结果 如 图 15 所 示 。 


EE 损伤 指标 刀 ， 
--- 预警 线 


起 0.6 
站 | ssecJ===============EE= 
志 0.4 
颈 
0.2 
0 
| 这 3 4 5 
梁 段 号 
(a) 实验 损伤 工 况 1 
3.0 
2.5| E 损 伤 指标 D， 
--- 预警 线 


(b) 实验 损伤 工 况 2 


图 15 实验 损伤 识别 结果 图 


Fig.15 Damage identification results of experiments 
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图 15(a) 中 只 有 2 号 梁 段 的 损伤 指标 明显 超过 
预警 线 ,对 应 于 实验 损伤 工 况 1 的 5 单元 损伤 ;实验 
损伤 工 况 2 为 多 损伤 工 况 。 图 15(p) 中 ,2 号 梁 段 和 
5 号 梁 段 的 指标 值 均 超过 了 预警 线 , 与 预 设 的 5 单 
元 和 11 单元 损伤 相符 ,实验 结果 与 预 设 工 况 吻 合 程 
度 较 好 。 实 验 结果 说 明 本 研究 所 提 方 法 具有 一 定 的 
实用 性 和 有 效 性 。 


a 


5 结 论 


本 研究 提出 了 一 种 在 移动 荷载 作用 下 基于 因子 
分 析 的 简 支 梁 桥 损伤 识别 方法 。 该 方法 采用 因子 分 
析 法 对 采集 的 位 移 响 应 进行 降 噪 处 理 ,并 结合 统计 
方法 和 数据 融合 方法 来 对 桥梁 结构 的 损伤 进行 定 
到 数值 分 析 和 模型 实验 均 验 证 了 该 方法 的 有 效 
得 到 以 下 结论 。 
1) 因子 分 析 可 以 将 采集 的 每 一 个 响应 数据 分 成 
公 莫 因子 和 特殊 因子 两 部 分 ,测量 噪声 属于 公共 因 
沪 部 分 ,通过 提取 响应 的 特殊 因子 可 以 有 效 降低 噪 
下 对 损伤 识别 产生 的 影响 。 

2 ) 研 究 了 噪声 、 车 辆 移动 速度 及 载 质 量 对 损伤 
误 唱 产生 的 影响 。 研 究 的 噪声 强度 范围 为 10% ~ 
30 色 ,速度 范围 为 5 ~ 33 m/s, 载 质量 范围 为 1 ~9 4， 
满足 实际 高 速 路 况 的 需求 。 研 究 发 现 ,噪声 和 载 质 
量 对 识别 结果 的 影响 较 小 ,而 速度 的 影响 较 大 。 通 
过 颖 制 干扰 率 曲线 发 现 速度 和 载 质量 的 影响 规律 相 
似 习 结合 不 同 桥梁 模型 参数 下 的 车 速 干扰 紊 曲线， 
归 睛 出 采用 本 研究 方法 可 以 准确 定位 结构 微小 损伤 
的 建议 速度 范围 为 8~22 m/s。 

”3) 本 研究 方法 只 需要 对 桥梁 结构 施加 移动 车 辆 
激励 ,并 且 适 用 的 速度 及 载 质 量 范围 较 大 ,便于 实际 
操作 。 对 于 单 处 损伤 和 多 处 损伤 ,本 方法 均 能 实现 
准确 定位 , 旦 具有 很 好 的 噪声 鲁 棒 性 ,有 一 定 的 实用 
价值 。 

本 研究 方法 依然 存在 一 定 的 局 限 性 , 因为 本 研 
究 不 管 是 数值 模拟 还 是 实验 验证 , 均 采 用 简 支 梁 桥 
的 模型 ,对 于 其 他 类 型 桥梁 的 损伤 识别 仍 待 进一步 
研究 。 此 外 ,本 研究 方法 只 能 识别 简 支 梁 桥 的 损伤 
位 置 ,从 损伤 识别 结果 也 可 以 看 出 ,本 研究 损伤 指标 
并 不 能 反映 结构 的 损伤 程度 ;所 以 ,下 一 步 还 需要 对 
结构 损伤 程度 的 量化 进行 研究 。 
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